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Abstract: We provide a unified description of fermion masses and mixing angles in the 
framework of a supersymmetric grand unified SO(10) model with anarchic Yukawa cou¬ 
plings of order unity. The space-time is five dimensional and the extra flat spatial dimension 
is compactified on the orbifold S^I{Z 2 x Z' 2 ), leading to Pati-Salam gauge symmetry on 
the boundary where Yukawa interactions are localised. The gauge symmetry breaking is 
completed by means of a rather economic scalar sector, avoiding the doublet-triplet split¬ 
ting problem. The matter fields live in the bulk and their massless modes get exponential 
profiles, which naturally explain the mass hierarchy of the different fermion generations. 
Quarks and leptons properties are naturally reproduced by a mechanism, first proposed 
by Kitano and Li, that lifts the SO(10) degeneracy of bulk masses in terms of a single 
parameter. The model provides a realistic pattern of fermion masses and mixing angles for 
large values of tan /3. It favours normally ordered neutrino mass spectrum with the lightest 
neutrino mass below 0.01 eV and no preference for leptonic CP violating phases. The 
right handed neutrino mass spectrum is very hierarchical and does not allow for thermal 
leptogenesis. We analyse several variants of the basic framework and find that the results 
concerning the fermion spectrum are remarkably stable. 
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1 Introduction 

The diversity of elementary particles and of their fundamental interactions observed at 
the energies probed in various experiments so far finds an elegant description in Grand 
Unified Theories (GUT). The unification of strong and electroweak interactions in GUT 
also leads to the unification of fundamental fermions. Such a unification can be partial or 
complete depending on the choice of unified gauge symmetry. As it is well known, one of 
the most attractive choice is the GUT based on the SO(IO) group [1]. All the Standard 
Model (SM) fermions of a given generation can be accommodated in the 16 dimensional 
spinorial representation of SO(IO), together with an additional fermion singlet under the 
SM gauge group. This new fermion can be identified as right handed (RH) neutrino, a 
partner of the weakly charged neutrinos in the seesaw mechanism of type I [2-6]. 

Like most of the interesting proposals of physics beyond the SM, GUT also suffer from 
drawbacks. The most serious of them is perhaps the fact that the GUT do not provide 
a unique way to get the observed diversity in low-energy physics from the unity imposed 
at high energy. In general the unified gauge symmetry can be broken down to the gauge 
group of the SM in several different ways. If a spontaneous breaking is realised, this 
requires the presence of scalar fields in large representations of the gauge group, allowing 
arbitrariness in the construction and also leading to problems like doublet-triplet (DT) 
splitting [7-9] and large enhancement of the unified coupling above the scale of grand 


- 1 - 




unification (Mqut)- Another difficulty arises in GUT due to the complete unification of 
the matter fields. The quarks and leptons exhibit different mixing patterns and it is not 
obvious how to reproduce this feature in a unified framework. The other aspect of the 
flavour puzzle is the hierarchy among fermion masses (see [10] for recent review on the 
status of flavour puzzle). In typical SO(IO) GUT [11-42] in 4D, the Yukawa couplings can 
vary in a huge range, 0(10“®) to 0(1), and no advantage is obtained over the SM in this 
context, as realized by several dedicated attempts of explaining the fermion mass spectrum 
in some simple SO(IO) models, see for examples [43-53]. Clearly, lots of improvements and 
efforts are needed to come up with a realistic and natural theory of flavour based on GUT. 

Some of the above issues can be addressed by implementing the program of grand 
unification in higher space-time dimensions [54-66]. First of all, by adding a new spatial 
dimension compactihed on an orbifold we can break the gauge symmetry by se¬ 

lecting appropriate parities of the gauge fields [54] . Only the gauge fields with even parity 
survive on the 4-dimensional fixed points (or branes) leaving the corresponding gauge sym¬ 
metry unbroken. In this way, the breaking of SO(IO) down to the Pati-Salam (PS) gauge 
symmetry [67], namely SU(4 )cxSU(2)lxSU(2)r, have been studied in [63-65]. Once the 
symmetry is broken through the boundary conditions, one has the freedom to introduce 
on the branes scalar multiplets transforming only under the unbroken symmetry. As it 
was shown in [54], this offers an elegant solution to the DT splitting problem. A second 
important aspect concerns the flavour problem of GUT, which can greatly benefit from 
the presence of an extra compact dimension. In the framework proposed by Kitano and Li 
in [60], an SO(IO) model in five flat space-time dimension (5D) is realised, with the extra 
dimension compactihed on jZ^- The three generations of matter helds are kept in the 
bulk and their bulk masses create exponential prohles for the corresponding zero-modes. 
The inter-generational mass hierarchies is explained by 0(1) fundamental parameters. The 
difference between the quarks and leptons is reproduced by spontaneous breaking of the 
SO(IO) symmetry into SU(5)xU(l)x through a bulk scalar multiplet. A complete and pre¬ 
dictive model based on this idea has been constructed in [66] , showing that fermion masses 
and mixing patterns can be successfully described in terms of fundamental parameters of 
0 ( 1 ). 

In this paper, we provide a merger of these two basic ideas. We construct a 5D SO(10) 
model with N=1 supersymmetry (SUSY) in which the extra dimension is compactihed on 
an orbifold /(Z 2 x Z 2 ) [68]. An N=1 SUSY in 5D is equivalent to N=2 SUSY in 4D [69]. 
The rehection under Z 2 breaks one of the SUSY while Z 2 is used to break SO(IO) down to 
the PS gauge symmetry. Thus the effective symmetry on one of the two branes is the PS one 
with N=1 SUSY. The further breaking of PS to the SM gauge symmetry is implemented 
by introducing appropriate helds on the brane. Fermions are described by 16 dimensional 
representations living in the bulk. As a consequence of the breaking of SO(IO) down to the 
PS symmetry the fermion zero modes fall into multiplets of the PS gauge group, namely 
(4, 2,1) or (4,1,2), depending on the Z 2 parity assignment, and a doubling of matter helds 
per each generation is required. This has the advantage of allowing different prohles for the 
zero modes of (4, 2,1) or (4,1, 2) in each generation. At this stage quark-lepton unihcation 
inherited from the PS symmetry still holds, and a new independent source of breaking of 
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the PS symmetry is required. This is obtained by the vacuum expectation value (VEV) 
of an adjoint scalar multiplet that spontaneously breaks SO(IO) into SU(5)xU(l)x giving 
rise to a distinct set of zero mode profiles. Such a breaking is flavour blind, introduces 
only one new parameter and contributes with different weights to lepton and quarks bulk 
masses. 

The model presented here provides a simple and viable alternative to the modified 
Kitano-Li (KL) model constructed by us in [66], based on the framework proposed in [60]. 
In comparison to that, the current model implements in a simpler way the GUT symmetry 
breaking and requires representations for the scalar fields with smaller dimensionality. The 
DT splitting problem does not arise since no color triplet is associated with the weak dou¬ 
blets introduced by us. The simplified scalar spectrum on the brane reduces the number of 
non-anarchic free parameters in the theory compared to the modified KL model, providing 
in principle a more predictive framework for the description of the fermion mass spectrum. 
While the number of independent parameters is still quite large, not allowing for precision 
tests of the model, we find that all fermion masses and mixing angles can be described with 
all the fundamental parameters of the theory of 0{\). A good agreement of the model with 
the data can only be obtained with large values of tan (3, where tan f3 is the ratio of the 
VEVs of two Higgs doublets used in the minimal supersymmetric standard model (MSSM). 
While both the normal and inverted ordering in the light neutrino masses can be obtained, 
the normal ordering is considerably less fine-tuned in the anarchic Yukawas. We derive pre¬ 
dictions for the CP violating phase in the lepton sector, the amplitude of the neutrinoless 
double beta decay and masses for the right-handed neutrinos. Within the same basic setup 
we also study another realisation of the Yukawa interactions, resulting in a model very 
similar to the modified KL model with an increased set of free parameters. A quantitative 
comparison of both the alternatives is also given. 

The organization of paper is as follows. We describe the model including the dynamics 
on bulk and on the branes in the next section. We then discuss how the fermion mass 
relations arise in the model in section 3. A qualitative comparison between the alternative 
models is given in this section. In section 4, we provide a detailed numerical analysis of the 
various options and discuss the results and predictions for the different observables. The 
study is finally concluded in section 5. 

2 An SO(10) model in five dimensions 

The model is based on a supersymmetric SO(IO) grand unified theory in five space-time 
dimensions [63-65]. The extra spatial dimension is compactified on an orbifold x ^ 2 ) 

where represents a circle of radius R. A periodic coordinate y parametrizes the circle 
and the action of the parity Z 2 {Z 2 ) is defined by y —)■ —y {y' —)• —y')-, where y' = y —7ri?/2. 
Points of the circle related by either Z 2 or Z 2 are identified. The interval between the two 
fixed points y = 0 and y = 7ri?/2 can be considered as the fundamental region. The other 
fixed points y = ttR and y = —t:RI2 are identified with the points y = 0 and y = 'kRI2, 
respectively. A generic bulk field y) can be categorized by its transformation properties 
under Z 2 x Denoting by P and P' the parities under Z 2 and Z 2 respectively, a field 
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(j)p^pi{x, y) with given parities (-P, P') can be expanded in terms of Fourier series as follows 
[65]: 


(t>++{x,y) 

(t>+-{x,y) 

(t>-+{x,y) 

4>—{x,y) 



( 2 . 1 ) 


Here n = 0,1, 2,... denotes the different 4D Kaluza-Klein (KK) modes of a given bulk field. 
In the free theory, upon the compactification, a 4D component (j)^{x) acquires a mass k/R, 
an integer multiple of the compactification scale 1/R. Only the field with (P, P') = (+, +) 

contains a massless mode and it is non-vanishing on both the branes. The field (j)-\ _ (4> _ 

vanishes on the y = 'kRI 2 (y = 0) brane, while (j) _vanishes on both the branes. 

The theory possesses N=1 SUSY in 5D which corresponds to N=2 SUSY in 4D [69]. 
We utilize the Z 2 symmetry to break N=2 SUSY down to the N=1 SUSY in 4D [68]. In our 
set-up, the matter and gauge helds propagate in the bulk. We introduce a 16-dimensional 
hypermultiplet 16-^ for each SM generation of fermions and 45-dimensional vector-multiplet 
45v under N=1 SUSY in 5D. In 4D, these correspond to a pair of N=1 chiral multiplets for 
16^ = (16,16'=), and a vector and chiral multiplets for 45v = (45^/, 45$). The breaking 
of N=2 SUSY down to the N=1 SUSY in 4D is achieved by assigning even Z 2 parity to 
the 16 and 45y multiplets and odd Z 2 parity to their superpartners 16'= and 45$. 

The Z 2 symmetry is used to break the SO(10) gauge symmetry down to the PS symme¬ 
try [63-65]. The PS gauge symmetry is isomorphic to SO(6)xSO(4) and hence the parity 
assignments with respect to P' should be appropriately chosen such that the generators of 
SO(6)xSO(4) remain unbroken. Under SO(6)xSO(4), the two index antisymmetric SO(IO) 
representation 45 decomposes as (15,1) -|- (1,6) -|- (6,4). The hrst two submultiplets are 
taken even and the last one is chosen odd under This assignment breaks SO(IO) down 
to the PS group and set to zero all the gauge fields, other than those of the PS group, 
on the y = ttR/2 brane. The gauge interactions on this brane respects only the PS gauge 
symmetry. On the y = 0 brane, the full 45y exists but only the PS gauge helds have 
massless modes. For these reasons, we call the y = TrR/2 brane “a PS brane” while the 
y = 0 brane “an SO(IO) brane”. 

Once the P' assignments for the gauge helds are chosen as above, the ones for the 
matter submultiplets follow from the invariance of the gauge interactions. Under the PS 
symmetry, the SO(IO) 16-plet decomposes as (4, 2,1) -|- (4,1, 2). It can be seen from the 
gauge interactions that (4, 2,1) and (4,1, 2) must have opposite P' charges. Therefore only 
one of the two possesses zero modes and is different from zero on the y = 7ri?/2 brane. To 
accommodate zero modes for a full SM fermion generation we have to double the 16-plet 
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5D N=1 

4D N=1 

4D N=1 in PS 

ip,n 


45^ 

45# 

(15,1,1) + (1,3,1) + (1,1,3) 

(+, +) 

45v 

(6,2,2) 

(15,1,1) + (1,3,1) + (1,1,3) 
(6,2,2) 

(+, “) 
(-,-) 

(“, +) 


16 

(4,2,1) 

(+, +) 

16 « 

16 '’ 

(4,1,2) 

(4,1,2) 

(+, “) 

(“, +) 


(4,2,1) 

(-,-) 


16 ' 

(4,2,1) 

(+, “) 


(4,1,2) 

(+, +) 

16 "’ 

(4,1,2) 

(-,-) 


(4,2,1) 

(“, +) 


Table 1. The parities P and P' of different SO(IO) multiplets and their Pati-Salam submultiplets. 


[63-65] and assign mutually opposite P' charges for the PS submultiplets. Therefore, we 
introduce 16]^^ per each generation in the bulk with P {P') equal (opposite) to that of the 
16 h. Notice that this doubling destroys the full quark-lepton unification achieved with 
only one copy of 16-plet per generation. We summarize the P and P' assignment of all 
the bulk fields in Table 1. 

We now discuss the symmetry breaking pattern in the model. The SO(10) sym¬ 
metry is broken down to the PS gauge symmetry on the branes by the action of 
We use the mechanism originally proposed by Kitano-Li in [60] to break the PS sym¬ 
metry down to the SU(3)cxSU(2)lxU(1)b-lxU(1)3r group. This can be achieved if 
an SU(5) singlet belonging to 45# develops a vacuum expectation value (VEV) which 
breaks SO(IO) into SU(5)xU(l)x in the bulk. The residual symmetry on the branes 
is SU(3)cxSU(2)lxU(1)b-lxU(1)3r which in turn has to be broken into the SM gauge 
symmetry by introducing appropriate 4D fields on the brane of interest. We will discuss 
the brane sector and the breaking of U( 1 )b-lxU( 1 ) 3 r down to U(1 )y later in this section. 
Let’s first discuss in details the dynamics in the bulk. 


2.1 The bulk 

The N=1 SUSY in 5D allows only gauge interactions in the bulk [68]. The 45# inter¬ 
acts with the chiral multiplets 16, 16', 16'’ and 16"’ through gauge interactions. The 
superpotential in the bulk is: 


Wbuik = 16 i 


m. 


i T dy '\/2(75 45# 


16i + 16''’ 


m[ + dy - V2g5 45# 


16' 


( 2 . 2 ) 


Here i = 1,2,3 denotes three generations of matter. The bulk masses can be chosen 
real and diagonal without loosing generality and are parametrized by rhj and m'. The 
invariance of Wbuik under Z 2 x Z 2 makes the bulk masses odd under both the parities 
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and they can be expressed as m = m sgn{y) and m' = m! sgn{y), where m and m' are 
real constants and sgn{y) has period irR. Performing a KK expansion for the matter 
fields, namely 16(x,y) = 16nix)fn{y), after the dimensional reduction one gets for the 
massless modes [60] : 

I 27T? 

fo{y) = Y e-'^y for 0 < y < 7ri?/2 . (2.3) 

The /o(y) is appropriately normalized in the interval [0,7ri?/2]. Similar expression for the 
profiles of the 16^ zero modes can be obtained by replacing m with m' in Eq. (2.3). The 
4D massless mode is localized at y = 0 (y = ttR/2) brane for positive (negative) value of m 
and its value is exponentially suppressed on the opposite brane. The exponential behaviour 
of the zero-mode wave-functions can be used to explain the hierarchies among the fermion 
generations. 

The bulk masses do not distinguish the profiles of quarks and leptons of a given gener¬ 
ation residing in the 16 or 16' and at this stage the observed differences in the quarks and 
lepton masses and mixing patterns cannot be reproduced. A very crucial correction to this 
picture can be achieved through the Kitano-Li mechanism [60]. The VEV of 45$ along the 
SU(5)xU(l)x direction introduces a correction to the bulk masses and distinguishes the 
profiles of the SU(5) submultiplets. As proposed in [60], this correction, which introduces 
a single new parameter, modifies the bulk masses according to 

ml = mi - \/2Qxg5{45c!>) , (2.4) 

where r = (10, 5,1) represent matter SU(5) representations and are the corresponding 
U(l)x charges: = —1, Qx — 3 and Qx = —5. The above modification in the bulk 

masses was argued to be able to generate viable hierarchies in quarks and leptons and this 
was demonstrated in a specific model [66] through a detailed numerical analysis. Expressing 
the dimensionful quantities in units of the cut-off scale of the theory A, we rewrite 

771 ^ 

< = ^ = /r, - Ql,kx , (2.5) 

where m = mi jA. and kx = -v/2y5(45$)/A. As discussed earlier, our Z 2 parity assignment 
allows massless modes for (4, 2,1) G 16, which contains the SM weak doublets of quarks 
and leptons {Q,L) and for (4,1,2) G 16' containing the weak singlet fields {u^, d'^, , N'^). 

The different matter fields within PS multiplets receive appropriate corrections from the 
VEV of 45$ proportional to their U(l)x charges: 

af = IJ'i + kx ] af = Hi - Zkx ; 

= y^'i + kx ] af = Hi - ^kx ; 

= y'i + kx ; af= Hi + ^kx ■ (2.6) 

In conclusion Hi and Hi responsible of splitting the prohles with respect to the PS 
submultiplets while kx with respect to SU(5) submultiplets. The zero mode prohles for 
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the various matter fields can be rewritten from Eq. (2.3) in terms of the dimensioless 
quantities as: 


nf{y) = \/A/o“ (?/) 



(2.7) 


where a = {Q, d^, L, e'^, N^) represents MSSM matter fields while c = AnR is a parame¬ 

ter which depends on the relative separation between the compactification scale and cut-off 
of the theory. 


2.2 The branes 

The N=1 SUSY in 5D forbids Yukawa interactions in the bulk which can be enabled on the 
branes by introducing a proper Higgs sector. As discussed earlier, on the y = TTRj2 brane 
only the PS gauge symmetry survives and one can introduce 4D fields filling representations 
of the PS gauge group. On the contrary, on the y = 0 brane full SO(IO) multiplets of 4D 
fields are required. Therefore the PS brane provides a more economical option in terms 
of the number of 4D fields. More interestingly, for light particles we can introduce only 
color singlet and electroweak doublet fields on the PS brane, avoiding the DT splitting 
problem. We introduce 4D chiral multiplets H , H ' transforming as (1,2,2), S ~ (4,1,2), 
S ~ (4,1,2) and T ~ (1,1,3) on the PS brane and 16^, IQh on the SO(IO) brane. The 
superpotential is 

Yi, 16 , 16 '^ + ^ + ... 

+ < 5 ( 2 /-^) u;^(i7,i7',S,S,r)-F5(y) u;o(16h,16//) , (2.8) 

where the first line in W corresponds to the Yukawa interactions responsible for the masses 
of matter fields, while Wt^ and wq are superpotentials for the chiral multiplets when the 
matter fields are turned off. The Y and Y' are complex 3x3 matrices while Yr is a complex 
symmetric matrix. Below we discuss the roles played by each of the brane fields. 

• S, S on y = 7ri2/2 brane 

These fields on y = 'i:R/2 brane play a multiple role. As discussed earlier, SO(IO) 
breaks down to SU(3)cxSU(2)lxU(1)b-lxU(1)3fi. One can construct two orthog¬ 
onal linear combinations of the generators of the two U(l)’s which can be identified 
with the generators of U(l)x and the SM hypercharge U(1)y- In our normalization 
convention, they read 

Qy = Y 3 R H-^— • (2.9) 

The fields S, S take a VEV along the U(1)y direction, trigger the breaking of 
U(1)b-lxU(1)3r down to U(1)y and contribute to the mechanism by which D-terms 
are canceled. The VEV of 45$ in the bulk generates D-terms on the branes [69-71] 


w = U.-f) i 
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associated to the U(l)x gauge symmetry. To preserve SUSY at high scale these D- 
terms have to be canceled by appropriate dynamics on the branes. The cancellation 
of the D-term on the y = TrR/2 brane can be achieved by the VEVs of S and S with 
the condition [69, 70]: 

= 2(45$) +55QI imf - |(S)|2) = 0 . (2.10) 

Here = — 5 is the charge under U(l)x of the component of S that acquires a VEV. 
Finally, the VEVs of S and S generate the masses for the right-handed neutrinos as 
shown in the first line in Eq. (2.8). 

• 16//, 16// on y = 0 brane 

The role of these fields on the y = 0 brane is similar to that of S and S on the other 
brane. The VEV of the singlet under SU(5)xU(l)x residing in 16//, 16// cancels 
the D-term on y = 0 brane if 


Do = -2(45$) TysQx (l(16//)P - \{ 16 h)\^) = 0 , (2.11) 

where = —5 is the U(l)x charge of the SM singlet in 16//. 

• H, H', T on y = ttR/I brane 

The H and H' are responsible for Dirac type masses of all the fermions. Each of the 
H and H' contains a pair of Higgs doublets which get mixed through the following 
terms in in Eq. (2.8): 

+ rnHH' + XTHH' + T{\hH^ + + ... (2.12) 

where dots stand for additional terms involving the S, S fields. Decomposing H and 
H' into electroweak doublets, H = (77^,77/) and 77' = (77(j,77^), one obtains the 
following mass term after the electroweak singlet in T acquires a VEV: 


(77, 77 ;) M 



with M = 


Mh m- X{T) 
m + X{T) Mjji 


(2.13) 


Here Mh,h' are redefined including the contributions coming from the VEV of T. 
All the mass parameters are assumed to be much heavier than the electroweak scale, 
possibly close to the GUT scale. One can arrange a pair of nearly massless Higgs 
doublets, by enforcing one eigenvalue of Ai being much smaller than the other. Such 
a pair would be an admixture of doublets residing in 77 and 77' and can be written 
as 

^u,d — COS T ^A\.9u,dH^ (i (^'1^) 

where, in the limit det(A7) = 0, the mixing angles read 


9u,d — 


1 _i 

- tan 

2 


/ 2MH'{mTX{T)) \ 
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The other combinations orthogonal to hu and obtain masses as large as the GUT 
scale. Below the GUT scale, the model contains only one pair hu,d which plays the 
role of MSSM Higgs doublets and triggers electroweak symmetry breaking. Glearly, 
getting much lighter than the GUT scale requires a fine-tuning of the parameters 
in (2.13). As we show in the next section, both H and H' with 9u 7 ^ Od are needed 
to generate viable quark mixing angles. Hence a non-vanishing (T) is required. We 
note that the VEV of T breaks SU(2 )r by keeping U(1)3 r unbroken and does not 
give any additional contribution to the D-terms on the PS brane. 

The model involves multiple scales of symmetry breaking. 

SO(IO) PS SU(3)cxSU(2)lxU(1)b-lxU(1)3r SM 

For simplicity, we take all these scales very close to each other and identify them with the 
GUT scale Mqut- Below the GUT scale the theory looks like the MSSM and we expect 
standard SUSY gauge coupling unification [72-74]. In order to suppress the higher order 
corrections in Eq. (2.8), we take c = AttR ~ 0(100) so that the cut-off of the theory, A 
can be lifted up to the Planck scale (see [66] for more discussions on the allowed range of 
the c parameter). The higher order corrections are at the percent level and remain smaller 
than experimental uncertainty in the fermion mass data we adopt. The theory provides a 
predictive framework for fermion masses and mixing angles, to be discussed in details in 
the following section. 

Before ending this section we notice that Yukawa interactions can also be present on 
the SO(IO) brane. A possibility is that all Yukawa interactions are localised at y = 0. 
In this case the dynamics on this brane becomes very similar to the one described in the 
modified Kitano-Li model discussed by us in [66]. The scalar content on the y = 0 brane 
in [66] consists of 10/^, 120/^, 126//, 126// and 45//. This combination of fields provides 
the most economical setup for viable fermion masses and mixing angles, a solution of the 
DT problem using the missing partner mechanism [8, 9, 75, 76] and a consistent GUT 
symmetry breaking. All these features are already discussed in details in [66] and we do 
not repeat them here. In the next sections we will briefly comment on the possibility to 
adopt the same scalar sector for the y = 0 brane in the present setup and we will study its 
potential in explaining the fermion masses and mixings. 

3 Fermion masses on the branes 

The bulk and brane superpotentials in Eqs. (2.2) and (2.8) encode the information about 
the fermion masses and mixing angles. As discussed earlier, the Z 2 parity and the VEV of 
45$ split the zero-mode profiles of various fermions, while the mixing of H and H' leads 
to the following effective 4D Yukawa couplings; 

yu = Fq Yu Fuc ; yd = Fq Yd Trfc ; y^ = Fl Y^ F^c and = Fl Yu F^c , (3.1) 

where yu,d,e,u stand for the 3x3 matrices of dimensionless Yukawa couplings of down-type 
quarks, up-type quarks, charged leptons and Dirac neutrinos, respectively. The profile 
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matrices are given by 


Fa 


(n^{-KR/2) 0 0 

0 rq{'KRl2) 0 

\ 0 0 nf(7ri2/2) 


with a = {Q, u‘^, F, L, e'^, N^) 


(3.2) 


where nf{y) are defined in Eq. (2.7). The arise from the mixing of MSSM-like Higgs 
doublets in H and H' and, from Eq. (2.14), can be explicitly represented in terms of 
fundamental Yukawas as follows: 


= cos Ou^dX - sin Ou^dX' . (3.3) 

Considering the fact that (Euc )33 ps 0(1) ^ (Fu<^) 22 , {Fu^)ii and the same for F^c, one 
obtains ~ FQYu^dY^ d^Q- ^ common Yukawa Yu = Yd leads to an unrealistic 

scenario of nearly vanishing quark mixing angles. Therefore we require (a) at least two 
pairs of Higgs doublets allowing for different Y and Y' and (b) unequal mixing 6u ^ Od 
to ensure that Yu and Y^ are different. The latter condition is satisfied in our model by a 
SU(2 )r triplet field T as shown in Eq. (2.14). After the electroweak symmetry breaking 
through the VEVs of hu^d, one obtains the mass matrices: 

^ V cos l3 Yd,e and Mu = v sin j3 Yu, (3.4) 

where tan/3 = (hu) / (hd) and v = a/ (hu)^ + (hd)"^ = 174 GeV. 

The RH neutrinos receive masses through the U(1)b-l breaking VEVs of S and are 
given as: 

F[r = vr Fatc Yr Fjsfc , (3.5) 

where vr = (S)^/A represents the seesaw scale. If the cut-off of the theory is raised to the 
Planck scale, the seesaw mechanism takes place two order of magnitude below the GUT 
scale, the right scale to generate viable neutrino masses. The light neutrinos gain masses 
through the type-I seesaw mechanism and their mass matrix can be expressed as 

Mu = (YuY^^Y^) Fl . (3.6) 

VR 

The model contains 24 complex parameters of 0(1) (9 each in Y and Y' and 6 in Yr) as 
the fundamental Yukawa couplings. In addition, it has two Higgs mixing angles 9u,d and 7 
bulk mass parameters y,i, and kx- 

As it was originally proposed in [60], the bulk masses in Eq. (2.6) can generate different 
hierarchies in Fq and Fr, which in turn explain the observed differences in the quark and 
lepton mixing patterns and mass hierarchies. The SO(IO) breaking by Z 2 distinguishes the 
profiles of left and right handed fields but it still maintains the quark-lepton unification. 
A milder hierarchy among neutrino masses and large lepton mixing angles result from the 
VEV of 45#, which distinguishes profiles of different SU(5) submultiplets within the 16 
and 16b This model differs from the one presented in [66] in the following ways: 
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• In comparison to [ 66 ], the current model has three more bulk masses. This provides 
more freedom in the profiles of zero-mode fermions. For example, the effective SU(5) 
symmetry in the profiles is broken once m* 7 ^ and, unlike in the previous model, 
one can distinguish between the masses of down-type quarks and charged leptons 
even if = Y^. 

• An important difference with respect to [ 66 ] is the simplification of the Higgs sector 
on the brane. In [ 66 ], consistent fermion masses and a solution of the DT splitting 
problem through the missing partner mechanism required lO/f -|- 120// Higgs repre¬ 
sentations, which contain three pairs of MSSM-like Higgs doublets. In the current 
model, only two pairs are required and this reduces the Higgs mixing parameters 
from eight to two. 

• The scalars introduced on the PS brane are in representations of smaller dimen¬ 
sionality compared to the brane sector fields in [ 66 ]. In particular, realistic Yukawa 
couplings only require a pair of (1, 2, 2) fields on the PS brane. The DT splitting is 
automatically solved since no colour triplets are present in the relevant Higgs multi- 
plets. However we need to arrange only one pair of light doublets and this requires 
an appropriate potential with a fine-tuning, as explained in the last section. 

As recalled at the end of the previous section, all Yukawa couplings can be also localised 
on the SO(IO) brane at y = 0. We can adopt the same scalar sector as in the model 
discussed in [ 66 ], remarking however a couple of differences with respect to our previous 
model. There are three more bulk masses in the current setup due to the doubling of matter 
fields in 16 and 16' and the Yukawa matrix Yio (Y 120 ) is not symmetric (anti-symmetric) 
in generation space, with several new parameters of 0{1). Clearly, this model does not 
provide any improvement in comparison to the old model as far as the field content and 
dynamics on the brane are concerned. It is however characterized by more parameters, 
which provide more flexibility in reproducing the correct pattern of fermion masses and 
mixing angles. We will provide a quantitative analysis of this improvement in the next 
section. 

4 Numerical analysis and results 

We now discuss in detail the viability of the model in explaining the observed data of 
fermion masses and mixing parameters and analyze its prediction for the observables which 
have not been measured yet. Our approach is similar to the one followed by us earlier in [ 66 ]. 
We take an idealized set of data for fermion masses and mixing parameters extrapolated at 
the GUT scale in the MSSM and check the viability of the model in reproducing them. As 
in [ 66 ], we use the results obtained in [77] for the charged fermion masses and quark mixing 
parameters. The extrapolation was carried out in the MSSM assuming a SUSY breaking 
scale of about 500 GeV, and for different values of tan/3. We perform the viability analysis 
for two representative values of tan/3, 10 and 50. After our previous analysis, the results 
of the global fit of neutrino oscillation data have been updated [78] taking into account the 
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Observables 

tan/3 = 10 

tan /3 = 50 

yt 

0.48 ±0.02 

0.51 ±0.03 

Vb 

0.051 ±0.002 

0.37 ±0.02 

Vr 

0.070 ±0.003 

0.51 ±0.04 

mulrric 

0.0027 ±0.0006 

0.0027 ± 0.0006 

md/rus 

0.051 ±0.007 

0.051 ± 0.007 

rrie/m^ 

0.0048 ± 0.0002 

0.0048 ± 0.0002 

mc/mt 

0.0025 ±0.0002 

0.0023 ± 0.0002 

mslrrib 

0.019 ±0.002 

0.016 ±0.002 

m^/m^ 

0.059 ±0.002 

0.050 ±0.002 

|K«| 

0.227 ±0.001 

lUf,! 

0.037 ±0.001 

\Vub\ 

0.0033 ± 0.0006 

Jcp 

0.000023 ± 0.000004 

As/10-5 eV^ 

7.50 ±0.19 

(NO or 10) 

A^/10-5 eV^ 

2.457 ±0.047 (NO) 

2.449 ± 0.048 (10) 

sin^ 012 

0.304 ±0.013 (NO or 10) 

sin^ 023 

0.452 ± 0.052 (NO) 

0.579 ±0.037 (10) 

sin^ 013 

0.0218 ±0.0010 (NO) 

0.0219 ±0.0011 (10) 


Table 2. The GUT scale values of the charged fermion masses and quark mixing parameters from 
[77] and neutrino masses and mixing parameters from an up-to-date global fit analysis [78]. NO 
(10) stands for the normal (inverted) ordering in the neutrino masses. 


most recent data available till the summer 2014. We take these updated low-energy values 
of neutrino mass squared differences and lepton mixing angles, neglecting RGE corrections. 
Such an approximation is valid if neutrino masses are hierarchical [79-81] and indeed this is 
realized in our model as we will show in this section. Following the widely adopted strategy 
in this kind of analysis [48-53], the data we use are the result of a specific extrapolation 
and should be taken as a representative set of GUT scale inputs. The actual data depends 
on features such as the SUSY breaking scale, SUSY scale threshold corrections, which can 
be estimated only when the exact mechanism of SUSY breaking is known [82-84]. Keeping 
these uncertainties in mind, we believe that if a given model can fit a representative set 
of data very well, then it will be able to reproduce with a similar accuracy and success 
the actual data, by slightly varying the underlying parameters. We summarize the various 
observables and their input values in Table 2. We employ minimization technique to fit 
the free parameters of the models with the data. See the details in [66] for the definition 
of the function and discussion on the optimization technique. 

4.1 Results for the PS brane 

We first analyze the Yukawa interactions on the PS brane. The compatibility of the model 
with anarchic Yukawa structure is tested in two ways. We first fit an idealized data set to 
the model by minimizing the x^ with respect to all the free parameters. The range of 0(1) 
Yukawa couplings is restricted to be \ Yij\, \Y-j\, G [0.5, 1.5] keeping the phases in the 
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full range [0, 27r]. The aim of this exercise is to assess whether our model can accommodate 
the data or not. We carry out this exercise assuming normal (NO) or inverted ordering 
(10) in the light neutrino masses and each of the two cases is analyzed for two values of 
tan (3. We get poor fits for small tan (3 corresponding to minimized values ~ 100 and 
~ 300 for NO and 10 cases respectively. The results for tan/3 = 50 are displayed in Table 
3 for which we get good fits for both NO and 10 cases. As it can be seen, all the data are 



Normal ordering 


Inverted ordering 

Observable 

Fitted value Pull 

Fitted value 

Pull 

yt 

0.51 

0 

0.52 

0.33 

Vb 

0.37 

0 

0.38 

0.50 

Vr 

0.51 

0 

0.51 

0 

mulrric 

0.0027 

0 

0.0028 

0.17 

rud/ms 

0.051 

0 

0.052 

0.14 


0.0048 

0 

0.0048 

0 

Wc/mt 

0.0023 

0 

0.0023 

0 

ms! mi, 

0.016 

0 

0.017 

0.50 

m^/rrir 

0.050 

0 

0.050 

0 

\Vus\ 

0.227 

0 

0.227 

0 

\Vcb\ 

0.037 

0 

0.037 

0 

\Vub\ 

0.0033 

0 

0.0030 

-0.50 

Jcp 

0.000023 

0 

0.000023 

0 

^s/^A 

0.0305 

0 

0.0305 

0 

sin^ 012 

0.304 

0 

0.304 

0 

sin^ 023 

0.452 

0 

0.442 

-0.20 

sin^ 013 

0.0218 

0 

0.0218 

-0.10 

Xmin 

« 0 

« 0.96 


Predicted value 


Predicted value 

[meV] 

3.9 


10.6 


\m0/d\ [meV] 

4.96 


48.2 


sin Sqp 

-0.39 


-0.89 


Mn, [GeV] 

190 


7.12 


Mat, [GeV] 

8.02 X 10^ 


6.75 X 10^ 


Mn, [GeV] 

1.43 X 1014 


1.38 X 1014 


VR [GeV] 

0.04 X lOi® 


0.056 X lOi® 



Table 3. Results from numerical fit corresponding to minimized x^ for normal (NO) and inverted 
ordering (10) in neutrino masses. The fit is carried out for the GUT scale extrapolated data given 
in Table 2 for tan /3 = 50. The input parameters are collected in the Appendix. 

fitted with negligible deviations from their central values. The model parameters obtained 
at the minimum of x^ are listed in the Appendix. The basic features of the best fit results 
are similar to the ones obtained in the previous model [66]. The observed hierarchies of 
quark and lepton masses requires \kx\ ~ l/^ 2 ,i| 5 1/^2 il ^ This in turn enforces a 

common bulk mass for quarks and leptons of the third generation and leads to approximate 
Yukawa unification yt ^ yt ^ yr, which prefers large tan/3 [85, 86]. 
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We now discuss the second kind of approach in which we do not fit the fundamental 
Yukawa couplings of the theory. We treat them as free 0(1) parameters and restrict their 
absolute values within the range 0.5 - 1.5, allowing arbitrary phases. For given values of 
these couplings, we minimize the function with respect to the bulk masses and Higgs 
mixing angles. We repeat this procedure many times, each time generating randomly a 
new set of Yukawa couplings. We fit 17 observables with respect to 9 free parameters (7 
bulk masses and 2 Higgs mixing angles), leaving v = 8 degrees of freedom (dof). The 
analysis is performed for tan/3 = 50 and for NO and 10 in the neutrino masses. The 
results are displayed in Fig. 1 where we plot the normalized distribution of the minimum 
jv. One can see a clear preference for the NO with respect to the 10. Even though one 
obtains a good best fit for 10 case in Table 3, this analysis shows that the solution requires 
more fine-tuning in the underlying Yukawas compared to the one obtained for NO. The 



X min/ 


Figure 1. The probability distributions of minimized jv for NO (blue) and 10 (red) in neutrino 
masses and for tan jS = 50. 

thresholds corresponding to a given probability value p and the number of cases satisfying 
the thresholds for different p—values are listed in Table 4. For p > 0.001, we find 0.5% 


p—value 

0.10 

0.05 

0.01 

0.001 

xLn (for 1^ = 8) 

< 13.36 

< 15.51 

<20.09 

< 26.12 

successful cases (NO) 

0.03% 

0.05% 

0.15% 

0.48% 

successful cases (10) 

< 10-3% 

< 10-3% 

< 10-3% 

0.005% 


Table 4. The rate of successful events obtained for different p—values from random samples of 
0(1) Yukawa couplings in case of normal and inverted ordering in the neutrino masses. 

cases providing the acceptable values of the Xmin — 26.12. The distributions of the bulk 
mass parameters and physical predictions for the NO case with p > 0.001 are given in Fig. 
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2 and 3 respectively. 
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Figure 2. The distributions of bulk mass parameters fitted with Xmin/^ < P > 0.001) in 

case of NO and tan^ = 50. The green (red) distribution corresponds to unprimed (primed) bulk 
mass parameters. 

One finds preference for positive bulk masses for the first and second generations, 
which are localized close to the y = 0 brane. The third generation is localized on the PS 
brane with a negative bulk mass. From the distributions shown in Fig. 2, it is clear that 
the SO(10) breaking by Z 2 , which distinguishes m and is crucial in generating realistic 
fermion masses in this model. This is particularly true for the first two generations where 
difference between yLi and /i) is significant. Notice that this difference is the only source 
of breaking of the mass degeneracy between the charged leptons and down-type quarks 
in this model. The kx parameter is required to be positive and of the order of the bulk 
masses of the first two generations. Among the observable quantities in the lepton sector, 
the lightest neutrino mass is predicted to be below 10 meV corresponding to strongly 
hierarchical neutrinos. The effective mass of the neutrinoless double beta decay lies 
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Figure 3. The Yukawa interactions of PS brane: prediction for various observables obtained for 
p > 0.001 (corresponding to Xmin/^ < 3.27 for i/ = 8) in case of normally ordered neutrino masses 
and tan/3 = 50. The black points in the bottom-right panel are model predictions while the green 
(red) regions are the allowed ranges for \mpp\ and the lightest neutrino mass in case of NO (10). 
The different horizontal and vertical grey bands correspond to the currently excluded regions by 
GERDA-I [87] and Planck Cosmic Microwave Background measurements and galaxy clustering 
information from the Baryon Oscillation Spectroscopic Survey [88]. The dashed lines indicate the 
near future reach of GERDA-II and KATRIN [89] experiments. 


in the range 1-5 meV, which is beyond the reach of the current generation of experiments. 
Future detection of neutrino masses well above 0.05 eV and/or of [m^/?] well above the 
range 1-5 meV would rule out the present model. Since the CP violation is coming from 
anarchic 0(1) Yukawas, we get no particular preference for the Dirac CP phase in the 
lepton sector. The model do not favour specific values also for the Majorana CP phases 
as revealed from the correlations between the and the lightest neutrino mass in the 

bottom-right panel in Fig. 3. 

Since the RH neutrinos are accommodated in 16-plets, their masses are predicted 
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once the masses and mixing angles of remaining fermions are fitted. The predictions are 
displayed in Fig. 4. The spectrum of RH neutrinos turns out to be very hierarchical. 



Figure 4. The Yukawa interactions of PS brane: prediction for the masses of RH neutrinos and 
vr = (S)^/A obtained for p > 0.001 (corresponding to Xmin/^ < 3.27 for = 8) in case of normally 
ordered neutrino masses and tan/3 = 50. 

This is a consequence of the large U(l)x charge of RH neutrinos which generates very 
large corrections in the bulk masses of the hrst and second generations making Ni ^2 niore 
sharply localized on y = 0 brane compared to the other fermions. Since kx ^ the 
third generation RH neutrino remains localized on the PS brane and one gets Mx^ ~ u/j = 
(E)^/A. We obtain relatively light spectrum for the first two generation RH neutrinos 
corresponding to Mx 2 £ 10^°] GeV and Mx-^ G [10^, 10®] GeV. This is in contrasts to 

generic 4D SO(IO) GUT models [51, 52] where they turn out to be relatively heavier. We 
also obtain the prediction for vp after correctly fixing the scale of solar and atmospheric 
neutrinos. This is shown in Fig. 4. One finds (S) ~ Mgut from the preferred values of 
vr which is of the same order as required by the cancellation of the D-term in Eq. (2.10). 
Note that |(S)| > |(E)| ~ Mgut is required since kx = \/ 2 fl' 5 ( 45 #)/A is positive. 

The spectrum of RH neutrinos is strongly hierarchical in our model. In the standard 
thermal leptogenesis [90] scenario, the final lepton asymmetry is dominated by the lepton 
number violating decays of the lightest RH neutrino. In this case the successful leptogenesis 
generically requires [91-93] 

Mjvi > 3 X 10® GeV . (4.1) 

Clearly, this condition is not respected in our model. To further assess the viability of this 
scenario, we perform a global fit imposing Eq. (4.1) in our model. We get Xmin ~ 
ruling out strongly the possibility of the Vi-dominated leptogenesis. An alternative is 
to consider N 2 or Vs-dominated leptogenesis, where the lepton flavour effects play an 
important role [94]. In this case, the lepton asymmetry is mainly generated by N 2 or 
V 3 decays. The lepton doublets produced in such decays get completely incoherent in 
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flavour space before the wash-out by the light RH neutrinos becomes active [95-99]. The 
wash-out acts individually on each flavour asymmetry and it is less efficient. In this case a 
certain combination of flavour asymmetry remains protected from the light RH neutrinos 
wash-out [94]. We have checked this possibility in our model using the best fit solution 
reported in Table 3 and in the Appendix. We find that N 2 is too light to create a sufficient 
asymmetry, while most of the asymmetry generated by is eventually washed out by 
N 2 and Ni, since these particles have sufficiently large couplings with lepton doublets and 
Higgs. Therefore, our preliminary investigations performed on the best fit solution indicate 
that leptogenesis cannot be successfully realized in this model. However a detailed analysis 
of this issue performing a global fit including the constraints imposed by flavoured 
leptogenesis would be required before ruling out leptogenesis in our model, which goes 
beyond the scope of the present work. 

4.2 Results for the SO(IO) brane 

We now investigate the naturalness of anarchic Yukawas on the SO(10) brane, as briefly 
discussed at the end of sections 2 and 3. The fermion mass relations are similar to the 
one already derived for the modified KL model in [66]. With respect to the modified KL 
model, we have three more bulk masses and several new Yukawa couplings in this model. 
We obtain good global fits for both NO and 10, when ion (3 = 50. Therefore we perform 
the second type of analysis in which we fit the 7 bulk mass parameters and 8 Higgs mixing 
parameters (see [66] for the details), by taking a flat random distribution for all the 0(1) 
anarchical parameters. The ranges of these parameters is chosen as in the previous case. 
Because of the new parameters coming from the Higgs mixing, with respect to the PS 
brane, we now have only v = 2 degrees of freedom. 

To compare this case to the previous one, we plot the distributions of for both 
of them and for NO in neutrino masses in Fig. 5. As it can be seen, both the distributions 
peak around similar values of jv- The SO(IO) case however has a relatively broader 
distribution leading to more successful cases for a given p-value. We get 7%, 15% and 30% 
successful cases for p-values greater than 0.05, 0.01 and 0.001 respectively (the correspond¬ 
ing thresholds for Xmin for zz = 2 dof are 5.99, 9.21 and 13.82). The substantial increase in 
the success rate in this case compared to that with Yukawas on the PS brane is attributed 
to the fact that we have six more mixing parameters providing more freedom in fitting 
the fermion masses and mixing angles starting from random Yukawa couplings. A similar 
improvements can be seen by comparing the success rates of this case with those of the 
modified KL model in [66]. The improved success rates in this case is due to three more 
bulk mass parameters, which allows better fitting of the data. 

The predictions for the various observables in the successful cases, corresponding to 
the p > 0.001, are displayed in Figs. 6 and 7. All the predictions are very similar to 
those obtained in the case of Yukawas on the PS brane and modified KL model in [66]. 
This shows that these predictions depend almost entirely on the dynamics of the bulk that, 
generating different zero-mode profiles, distinguishes the various fermion sectors. On the 
contrary, details of the brane interactions affects only very mildly our results. The main 
difference arising from the brane interactions in the different cases is the number of free 
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0.20 



Figure 5. A comparison between the Yukawa interactions on PS {y = TTRj2) and SO(IO) {y = 0) 
branes. The distributions are obtained for the normal ordering in the neutrino masses and for 
tan P = 50. 


0(1) parameters and Higgs mixing parameters. Our study shows that when the number 
of bulk mass parameters and Higgs mixing parameters increases also the rate of success, 
normalized to the number of degrees of freedom, increases. 

5 Conclusion and discussion 

Grand unified theories, proposed more than forty years ago, provide an elegant synthesis 
of electroweak and strong interactions, which greatly clarifies some of the crucial aspects 
of the SM such as particle classification, gauge anomaly cancellation, quantization of the 
electric charge and diversification of the gauge coupling constants. In SO(10) grand unified 
theories one fermion generation fits in a single representation of the gauge group, leaving 
room for a right-handed neutrino, which naturally gives rise to neutrino masses through 
the see-saw mechanism. This impressive feature is at the heart of the well-known problem 
of finding an acceptable description of quark and lepton masses and mixing angles, which, 
in the low-energy data, do not reflect at all such a complete particle unification. While it is 
certainly possible to accommodate the observed fermion spectrum by exploiting the most 
general Yukawa interactions allowed by the theory, not much is gained with respect to the 
SM since a huge hierarchy in the Yukawa coupling is needed to reproduce the data. 

An attractive framework where all the fundamental Yukawa couplings are of order 
one can be realized, even in SO(IO) grand unified theories, through the localization of the 
profiles for the zero-mode fermions in an extra dimension. Yukawa interactions are defined 
on one brane and the hierarchy among fermion masses of different generations depends 
exponentially on the bulk fermion masses. Quarks and leptons can be further differentiated 
by inducing a breaking of the SO(IO) symmetry in the bulk mass parameters. While the 
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Figure 6. The Yukawa interactions on SO(IO) brane: prediction for various observables obtained 
for the successful cases corresponding to p > 0.001 (or Xmin/’^ < 6-91 i/ = 2) in case of normally 
ordered neutrino masses and tan/3 = 50. See Fig. 3 for detailed description. 


generic ingredients of this construction are well-defined, a considerable freedom is left 
in model building, depending on the specihc implementation of the idea. In a previous 
work we relied on a spontaneous breaking of the grand unified symmetry, at the cost of 
introducing large SO(IO) representations for the symmetry breaking sector with a non¬ 
trivial mechanism to solve the doublet-triplet splitting problem. In the present work we 
have fully exploited the capabilities of the higher-dimensional construction, which allows 
for gauge symmetry breaking through compactification and offers a more economic solution 
to the doublet-triplet splitting problem. Since, compared to our previous model, the new 
construction significantly alters the allowed bulk masses and the Yukawa interactions, we 
think it deserves an accurate study of its properties, to assess whether the description of 
fermion masses and mixing angles remains the same or it undergoes major modifications. 

We propose a supersymmetric SO(10) model formulated in five dimension. The extra 
dimension is compactified on an orbifold /(Z 2 x Z 2 ) and plays a key role in break- 
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Figure 7. The Yukawa interactions on SO(IO) brane: prediction for for the masses of RH neutrinos 
and vr = (126//) obtained for the successful cases corresponding to p > 0.001 (or X^in/^ < 6.91 
for z/ = 2) in case of normally ordered neutrino masses and tan/3 = 50. 


ing the symmetries of the model. The compactification breaks N=2 SUSY down to 
N=1 SUSY in 4D and, at the same time, breaks SO(IO) down to the Pati-Salam group 
SU(4)c xSU(2)lxSU(2)r. A further reduction of the gauge symmetry is realized sponta¬ 
neously, through a symmetry breaking sector including an SO(IO) adjoint, automatically 
present in this 5D construction, and additional brane multiplets included with the purpose 
of canceling the D-terms of the theory. Below the GUT scale the residual gauge symme¬ 
try is that of the SM, which can be finally broken down to SU(3)cxU(l)em by a set of 
electroweak doublets localized on the PS brane. Matter multiplets, introduced in 16 rep¬ 
resentations of the GUT group as bulk fields, develop prohles for the zero-modes that are 
localized in specific regions of the extra dimensions. A different localization for the zero¬ 
mode profiles of the SM fermions is achieved by different bulk masses. As in the original 
Kitano-Li model, a universal parameter, proportional to the VEV of the adjoint of SO(IO), 
allows to distinguish the different SU(5) components inside a 16 representation. Moreover 
our framework allows for independent bulk masses for electroweak singlets and doublets of 
the various generations. Yukawa interactions can be localized either on the SO(IO) or on 
the PS brane. While we briefly commented on the first possibility, in our study we mainly 
concentrated on the PS case, since it offers the possibility of introducing an economic Higgs 
sector, which in particular automatically solves the DT splitting problem. 

Our model, with Yukawa interactions on the PS branes, has seven parameters con¬ 
trolling the bulk masses and two Higgs-mixing parameters, plus a large number of 0(1) 
Yukawa couplings. By fitting an idealized set of data, extrapolated at the GUT scale from 
the observed fermion masses and mixing angles, we find that the agreement is not trivial 
and requires a large value of tan/3. Moreover the case of inverted ordering in the neutrino 
mass spectrum requires much more fine-tuning in the Yukawa couplings than the case of 
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normal ordering. The lightest neutrino mass is predicted to be below 10 meV and the 
effective mass of the neutrinoless double beta decay lies in the range 1-5 meV. The 

model can be falsihed by the observation of either a non vanishing neutrino mass at KA- 
TRIN [89] or at the next generation of experiments. We hnd no preference for the 

Dirac CP phase of the lepton sector and the spectrum of RH neutrinos is predicted to be 
very hierarchical, which unfortunately is incompatible with the generation of the observed 
baryon asymmetry through thermal leptogenesis. 

It is remarkable that all these features remain essentially unchanged in several ver¬ 
sions of the SUSY SO(IO) model in 5D, having in common the property of describing the 
fermion spectrum through a set of zero-mode profiles able to distinguish the three gener¬ 
ations and the different SU(5) components inside a 16 representation. All the remaining 
features of the model such as the number of independent Yukawa couplings on the branes, 
the number of Higgs mixing parameters, the additional possibility of distinguishing weak 
doublets and singlets through the bulk masses, seem to play a secondary role which, at 
most, can influence the success rate of the model when statistical tests are performed. We 
conclude that the results are rather robust against modifications of the basic framework. 
We hnd very interesting the possibility of combining in a realistic scheme the anarchy of 
the underlying Yukawa couplings with the unihcation of one fermion generation implied by 
the SO(IO) GUT. On the weak side, as all models based on a large number of independent 
0(1) parameters, it is not possible to plan precision tests of these ideas to fully exploit the 
accuracy of existing data. 
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A Parameters obtained for the best fit solutions 

We provide the set of input parameters obtained for the best ht solutions corresponding to 
normal and inverted neutrino mass spectrum and tan/3 = 50 as presented in the Table 3. 


A.l Normal ordering 

The values of the Yukawa matrices and bulk masses appearing in Eqs. (3.1, 3.4, 3.6) at 
Xmin ~ 0 are as the following. We have removed some unphysical phases by redefining the 
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fields. 




Yd 


Yr 


/ 0.55863 0.94275 1.23911 

0.74927 1.49374 0.66883 

0.50804 e°'22i3ii 0.50000 1.26156 

/ 0.64691 0.71998 e-°-®i349i 0.52244 e2-'^284H \ 

0.80610 0.57351 e°-23467j 0.50398 

^1.01632 0.59252 0.63639 6-2-92490*^ 

/ 1.10716 0.70519 e0-94555i 0.81595 e-0'^'527ii \ 

0.70519e°-^'^®^^* 1.30773 e2-93543i i_o7719 

(^0.81595e"°-^®27ii q 77 ;L 9 g-o.i74iii 0.71443 y 



The corresponding bulk mass parameters are: 


{/ii, /r 2 , /Us} = {0.049590, 0.020895, -0.139245} , 
{/r'l, fi 2 , ^i 3 } = {0.066244, -0.013373, -0.463361} , 
kx = 0.042394 . 


(AT) 


(A.2) 


From the above parameters the profile matrices in Eq. (3.2) for various SM fermions 
can be expressed in terms of powers of the Cabibbo angle A as below. 


Fq = X 


Fl = X 



F,ip — Fpc—X 



- W 


Fx<^ = X 


(A.3) 


A.2 Inverted ordering 


The values of the Yukawa matrices and bulk masses appearing in Eqs. (3.1, 3.4, 3.6) at 
Amin ~ 0.96 are as the following. We have removed some unphysical phases by redefining 


the fields. 


t; 


Yd 


Yr 


/ 1.05063 e"2-27438i 0.50197 0.50108 

1.28888 0.95572 0.95749 

y 1.32079 ei-96363i 0.84379 1.03615 

/ 0.51388 e-2-467i9i 0.50192 0.72278 ei'99274i \ 

1.47850 e-^-34548i 0.63988 0.62270 

y 0.68440 e"^-^2037i 0.52781 e^-®2283* 0.506 1 8 e^-°3^28i j 

( 1.32057 1,34754 g-2.56275i 0.62345 e^-^2638i \ 

1.34754 e-2-56275i 1,44530 ei-87202i 0.57696 
0.62345 ei-^2638* 0.57 6 96 0.62830 e2-3ii8i* j 



(A.4) 
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The corresponding bulk mass parameters are: 


{/ri, // 2 , ^ 3 } = {0.056934, 0.023583, -0.212866} , 

{/i'l, fj.2, /U 3 } = (0.088673, -0.025229, -0.421995} , (A.5) 

kx = 0.045419 . 


From the above parameters the profile matrices in Eq. (3.2) for various SM fermions 
can be expressed in terms of powers of the Cabibbo angle A as below. 


Fq = ^ 


Fl = X' 



Fm = 


Fxc = a ' 


The — FLc—A 


0.1 



(A. 6 ) 
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